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CN " Abstract 

In this article we show and implement a simple and efficient method to strictly 
locate eigenvectors and eigenvalues of a given matrix, based on the modified cone 
condition. As a consequence we can also effectively localize zeros of complex poly- 
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cj ! 1 Introduction 

1.1 Motivation 

| Determination of the eigenvalues and eigenvectors of a matrix is important in many 

areas of science, for example in web ranking [8, 15], computer graphics and visualization 
[24], quantum mechanics [5], statistics [12, 20], medicine, communications, construction 
vibration analysis [3, 23]. 

One of the most important numerical methods designed to calculate the eigenvalues 
and eigenvectors of matrix A is the power method [18, 21]. It is used to determine a 
maximum module eigenvalue of A and the corresponding eigenvector v. The limit of 
the product , where w is a randomly chosen element, is the vector corresponding 

to the largest eigenvalue. The eigenvalue can be calculated from the Rayleigh quotient 
A = ^j-^y- The most common method of solving the full eigenvalue problem, i.e. finding 
all the eigenvalues and corresponding eigenvectors, is the QR method [6, 13]. 

There are methods for locating eigenvalues such as Gerschgorin theorem [7] from 
1931. This theorem allows to strictly locate the position of the eigenvalues of the matrix 
with real or complex coefficients. However, it does not allow to localize the eigenvectors. 

With the growing importance of these concepts there is a need to look for new 
methods of localizing simultaneously eigenvalues and eigenvectors. We do not know 
strict and efficient methods for locating eigenvectors of real or complex matrices, so we 
want to fill this gap. Our aim is to create and analyze a new method of strict location 
eigenvectors and eigenvalues with the use of interval arithmetic 1 . Using the fact that 
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x For language CH — h one can use libraries such as 'boost' [1] or 'CAPD' [3]. 
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the polynomial W(x) = x n + a n -\x n 1 + . . . + aix + ao of the n-th degree equals the 
determinant of the matrix 

" 1 ... 

1'-. : 

: : '•• 1 
... 1 

_— ao —at — a2 . . . — a n _i_ 

we show that we can effectively localize zeros of real or complex polynomials (for non- 
strict methods see [10, 16]). 

The content of this paper can be described as follows. In the following short subsec- 
tion we present an algorithm which is the final outcome of our theory. In the Section 
2 we introduce notation and present the simple properties of our theory, which is based 
the concept of cone condition. The notion of cone condition originally appeared in the 
late 60's in the works of Alekseev, Anosov, Moser and Sinai [22]. These techniques are 
used in the examination of differential equations [4, 9, 14, 19, 25]. In Section 3 we show 
the basic properties of operations. In the following section we present the main result 
of this article Corollary 4.1, which allow as a consequence obtain strictly localize the 
position of the eigenvector corresponding to the given eigenvalue in a small cone. In the 
last section we present a simple algorithm and some numerical applications of the theory 
on a few basic examples. In particular for a random matrix of size 5x5 with randomly 
generated values from the set [—1,1] we obtain the eigenvectors and eigenvalues with 
typical accuracy e ~ 10 -11 (see Example 5.1) 2 . 

1.2 Basic Algorithm 

In this subsection we present a basic idea how our method can be used (more detailed 
explanation is given in last section). This algorithm determines the accuracy for nu- 
merical approximation of the eigenvectors and the corresponding eigenvalues. Now we 
introduce the following notation for clarity and simplicity of the presentation. 
A matrix A of dimensions N x N with elements aij E K is denoted by 

A := Kj]i< M -<jv G Mjv(K). 

For K,J C {1, ... ,N} by Awn := [akj]keK,jeJ we denote the sub-matrix with rows 
k £ K and columns j G J. We put m : n = {m, . . . , n} for m,n £ {1, . . . , N}. For 
k € {1, . . . , N} by ^ k we understand the set {1, . . . , N} \ {k}. The identity matrix is 
denoted by /. 

By [x] we understand interval representation of x G R that is an interval with rep- 
resentable ends such that x G [x]. For complex number z £ C we put [zj := [Re(z)] + 
[Im(z)]i, and for a matrix A = [a,ij]i<ij<N, {A} denotes the matrix [[<%]] i<i,j<jv- We 
refer the reader to [17] for more information concerning interval arithmetic. 

As a sup norm ||-A||oo we take the maximum sum of absolute values of the elements 
in rows of matrix A, that is 

N 

1 1^4 1 |oo = ||[ a ii']||oo := max y 1 1 ay | |oo- 

Ki<N *■ — * 
" " j=l 

2 In our program we use variables of type double. One could decrease the error by using numbers of 
arbitrary precision. 
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We are given a matrix A = [a^ ]i<ij<N & Mat(C) with pairwise disjoint single 
eigenvalues, Aj 7^ Xj for i ^ j and j,j g{1,..., iV}. We assume that we have calculated 
the numerical approximation of eigenvectors xi, . . . ,xn and eigenvalues 3 Ai, . . . , Ajv of 
the matrix A. To locate the exact position of A;-th eigenvector x k of A one can use the 
following algorithm. 

Algorithm 1.1. 1) Put P := [xi, . . . ,x N ], 

2) calculate interval hull [P _1 ]] of inverse of P (we call it the inverse interval matrix 
ofP) 4 , 

3) calculate the interval matrix J = {P' 1 } ■ {A} ■ [P], 

4 ) seek a constant e > as small as possible so that the matrix 
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(1.1) 

where the number 1 is at the intersection of the k-th row and k-th column, satisfies the 
following condition 



\Bf, 



- X 11 — x 

1 1 00 



\ B &k,k]\\o°- 



N 



.. } [k,l:N]\\oc < 

Then the fc-th eigenvector of A satisfies 

N 

x k G [x k ] + {-£,£} ■ ^[Re(xO] + l-e,ej ■ ^pm(xi)]i. 

We estimate the eigenvalue A^ from the Rayleigh coefficient 

^ € (^4xk,x k ) 
(xk,x k ) 

We show more precise details in the last section of this paper. 

Remark 1.1. Our method can be modified to work with multiple single eigenvalues 
(see Proposition 2.1). However, it does not work with multiple eigenvalues of the form 
(1.2) but this is a typical problem of most algorithms calculating the eigenvalues and 
eigenvectors. 



(1.2) 



3 To calculate the approximation of eigenvectors and eigenvalues one can use any numerical method. 
4 By [P^ 1 ] we denote such an interval matrix that Q _1 G [-P -1 ] for every Q e \P\. 
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2 Cones and Dominating Maps 

In this section we modify the concept of cones from [14], which will allow us to locate 
the eigenvectors of the matrix. 

Definition 2.1. Let E be a finite-dimensional Banach space with semi-norms ) • ( (we 
call it contracting), (•) (which we call expanding) such that 

|||x||| := max( )x( , (x)) 

defines an equivalent norm on E. By the r-norm for r > on the cone-space E we take 

|||x||| r := max( )x( , r • (x)). 
Definition 2.2. Let E be a cone-space. We define the r-contracting cone in E by 

)E{ r := {x G E : >x( > r(x)}, 

and the r -expanding cone in E by 

(E) r := {x G E : >x( < r(x)}. 

Note that 

E=)E( r U(E) r . (2.1) 

In the same way we define r-contracting cone and r-expanding cone in subspace E. If 
r = 1 we will omit the subscript r, in particular we speak of contracting cone. We 
introduce the scaling by r of semi-norms to have better control over size of the cones 
(see Figure 1), and which consequently will allow us to better locate the eigenvectors. 

(a) The contracting cone inlxl 




(c) The expanding cone in R x R. (d) The 1/2-expanding cone in R x R. 

Figure 1: The cones in the cone-space IxK. 

For a product E = E\ x E 2 we introduce cone-space (E, ) • ( , (•)) where we put 

)x( := ||xi||, (x) := ||x 2 || for x= (xi,x 2 ) G E\ X E 2 . 

Analogically, we define cone-space for direct sum E = E\ © Ei. 

In our main result, Corollary 4.1, the following proposition will play a crucial role. 




(b) The 2-contracting cone in R x R. 
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Proposition 2.1. Let E = E\ x E2 be cone-space such that dimFi = n, cUiilE^ = m 

and let r > be given. Assume that we have direct sum decomposition E = V± © Vi such 
that 

ViC)E{ r and V 2 C (E) r . 

Then 

dimFi = n and dimV2 = m. 

Proof. First we show that dirndl < n. For an indirect proof, assume that dimVi > n. 
Then there exist linearly independent vectors vi, . . . , v n+ i G V±. Obviously Vj = (wi, z\) 
for i G {l,...,n + l} and uniques Wi G E\, Zi G E<i- Since wi, . . . , w n+ \ G E\ and 
dimE 1 ! = n there exist a set of n + 1 scalars, ai, . . . , a„+i, not all zero, such that 

aiwi + . . . + a n+ iw n+ i = 0. 

Note that 

z := a\z\ + . . . + a n+ iz n+1 / 0, 

because otherwise vectors vi, . . . , v n+ i would not be linearly independent. Consequently 
we obtain 

(n+l n+1 \ 

i=l i=l / 

and thus r||z|| < ||0||, which implicate z = 0. We get a contradiction with the fact the 
sequence of vectors vi, . . . , v n+ i is linearly independent. 

The proof that dimV2 < m is analogous. Finally since dimE 1 = n+m and dimVi < n, 
dimV2 < m we obtain 

dimlq = n, and dimV2 = m. 

□ 

By an operator we mean a linear mapping between cone-spaces E and F. We denote 
the space of all operators by C(E, F). If F = E, we denote C(E, E) by C(E). 

Let A G C{E,F). We define 

)A( r := inf{i? G M + | ||| Ax\\\ r ^R\\\x \\\ r for all x G E : Ax G )F( r }, (2.2) 
(A) r := sup{i? G M+ I HI Ax\\\ r ^ R\\\x \\\ r for all x G E : x G (2.3) 

Let E C E, F C F be subspaces and let A G C(E, F) such that A{E) C F. We 
define 

)A\ z( r := inf{i? G M+ | ||| Ax||| r ^ R \\\ x ||| r for all x G E : Ax G )F( r }, 
(A\g) r := sup{i? G M+ I HI Ax\\\ r ^ i? ||| x ||| r for all x G E : x G (E) r }. 
Definition 2.3. We say that A G C(E,F) is r-dominating, if 

)A( r < (A) r . 

By V r (E,F) we denote the set of all A G C(E,F) which are r-dominating. If F = E, 
we denote the space V r (E,E) by V r (E). 

Observation 2.1. Let E C E, F C F be subspaces and let A G C(E,F) such that 
A(E) C F. We have 

)A\ g( r < )A( r and (A) r < (A\ & ) r . 
Moreover, if A G T> r (E, F) then A G £> r (£, F). 
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Proof. It is consequence of (2.2), (2.3) and Definition 2.3. □ 
Theorem 2.1. Let A G T> r (E,F) and let v G E be arbitrary. Then 

v G (£) r ylv G (F) r , 

Av G }F( r =4v£ )£( r . 

Proof. The proof is a simple modification of the proof of [14, Proposition 2.1]. □ 

Proposition 2.2. Let A G V r (F, G) and B G T> r (E, F). Then AoB G £> r (£, G) and 

)A o B( r < )A( r ■ )B{ r , (A o B) r > {A) r ■ {B) r . (2.4) 

Proof. To prove the first inequality from (2.4), let x G E and A o -B(x) G )G( r . From 
(2.2) and Theorem 2.1 we know that Bx G )F{ r . We have 



Hence 



\AoB(x)\\\ r <)A( \\\Bx\\\ r <)A( )B{ |||x 



)AoB(<)A{.)B{. 



Using (2.3) and Theorem 2.1, we obtain the second inequality from (2.4). 
As a simple consequence of (2.4) we obtain A o B G T> r (E, G). 



□ 



3 Operator Norms 

In this section we show how given an operator A to estimate )A{ r , {A) r . 

Consider two cone-spaces E = E\ X Ei and F = F\ x F2. Let A: E — > F be an 
operator given in the matrix form by 



A 



By 



A21 A22 

A\\\ r := max (pull + -||Ai 2 ||, r||A 21 || + ||A 22 | 



(3.1) 



we define the r-norm of operator A. Observe that it satisfies 

|||^4x|||. r < \\\A \\\ r ■ III x ||| r for x G E. 

Note that in general it is not, if E\ is not one dimensional, the classical operator norm 
for HI • || r . 

Observation 3.1. Let r G (0,oo), A G C(E 1 x E 2 ,F 1 x F 2 ) We put 



A 



An A12 
A21 A22 



and R 



L 
rL 



Then 



\A\ 



\Rx HI for x£fi, 
\RAR- 1 III . 



Theorem 3.1. Let A = [^]i<ij<2 G C r (Ei x E 2 ,Fi x F 2 ). 
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i) We have 

)A( r < ||Aii|| + -||Ai 2 ||. 

r 

ii) Additionally, if A22 is invertible, then 

(A) r > WA^ 1 ^ 1 - r\\A 21 \\. 

Proof. We show the proof for r = 1. For arbitrary r £ (0, 00) the proof needs a simple 
modification, see Observation 3.1. 

For the proof of the first inequality, we take x = (xi,x 2 ) £ E\ x E 2 such that 
Ax £ )F( . From Definition 2.2 we have 

+ A12X2II > ||^2i^i + -A22Z2II, (3.2) 

therefore 

lll^xlH = max(PnXi + A12X2II, \\A21X1 + A 22 a; 2 ||) 

( = 2) \\A n xi +A12X2II < {pull • ||xi|| + p 12 || • ||x 2 ||} 
< (Pn|| + P 12 ||)-|||x |||. 

For the proof of the second inequality, suppose that x = (xi,x 2 ) £ (E), where 
Xi £ Ei, x 2 £ E 2 . Then 

ll^ill < \\ x 2\\ = ||| x III • (3.3) 

We know that 

P22X2II > P^II^IHI > 0. (3.4) 

Finally, we obtain 

|||Ax||| > \\A 2l xi + A 22 x 2 || > P 22 x 2 || - P21X1II 

(3.4) (3.3) 

> ii^h^ihi- iia 2 i||||xi|| > (p^ir 1 - P21I1) • iiixiii . 



□ 



Example 3.1. Let us verify that the matrix A £ £(R x 

A = 



2 1.5 
1 5 



is dominating. 

By Theorem 3.1 we have 

)A( < 3.5 < 4 < (A), 

and therefore A is dominating. 

4 The Main Results 

In this section we show the main results of the paper, concerns the strict location of the 
eigenspace. 

We denote the spectrum of the operator A by cr(A) := {A £ C : A — XI is singular}. 
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Lemma 4.1. Let A G P r (P). T/ien 

A £ (7(A) =► |A| C [0, )A( r ] U [(A) r ,oo). (4.1) 
Moreover [0, )A( r ] n [(A) r , oo) = 0. 

Proof. Since A G V r {E) we get [0, )A( r ] n [(A) r ,oo) = 0. 

Now we show implication (4.1). Let A be an eigenvalue of A and let x G P be a 
corresponding eigenvector. By (2.1) we know that x G )P( r U(P) r . We consider two 
cases. First suppose that x G )P( r . Since x is an eigenvector, Ax = Ax, and thus 
Ax G )P( r . By (2.2) we get 

|A|<)A( r . 

Now suppose that x G (E) r . By (2.3) we get 

|A| > (A) r , 

which complete the proof. □ 

Let E be a vector space over the field K and let operator A: E — > E be given. One 
can easy deduce from the theorem of Jordan, see [11, Appendix to Chapter 4] for the 
general case, that if a = o~\ U o<i (for IK = R we assume additional that o\ = o~\ and 
&2 = 02) then there is a unique direct sum decomposition E = E ai © E U2 such that 
A(P CT1 ) C E ai , A(P CT2 ) C E a2 and a{A\ Eai ) = o"i, <r(A\ E „ 2 ) = a 2 . For c> we define 

P< c := P{ A . | A |< C } and P> c := P{ A . |a|> c }- 

Theorem 4.1. Lei P be a finite- dimensional cone-space and let A G T> r (E). Then there 
is a direct sum decomposition 

P = P<)A( r ® E >{A) T an d ^<>A( r C)P( r , P>(A) r C (P) r . 

Proof. For the clarity of presentation we consider only the case r = 1 and we omit the 
subscript r. 

From Lemma 4.1 and Jordan Theorem we obtain a direct sum decomposition 

E = P<)A{ © E>{A)i 

such that 

°( a \e< )M ) = [0, >A(] and a(A|^ (A) ) = [(A), 00). 

Now we show P< \m C )P( . Consider an arbitrary x G P< V4{ . The case when x = 
is obvious. Assume that 0. Without loss of generality we can assume that ||x|| = 1. 
For an indirect proof, assume that x ^ )P( . Then x G (P) (by 2.1). Let e > be 
arbitrary. From the fact that x G P< \m , we know that 

limsup v/|||A n x||| <)A( , 

n— >+oo 

and thus there exists m G N such that for all n G N 

n > m => y\\\A n x\\\ < )A( +e. 
Since x G (P) and from Theorem 2.1 we obtain 

x G (P) =>• Ax G (P) => • • • => A n x G (P). 
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Using (2.3) we get 

IPx|||>(A> |||x|||, 
|P 2 x||| = \\\A(Ax)\\\ > (A) |Mx|||>(A) 2 |||x| 

|P"x||| > {A) n HI x HI . 

Finally we have 



(A) = < V¥^xJ\ < )A( +e. 

Since e was arbitrary, we get a contradiction with the fact that A is r-dominating. 

Analogously to the proof of the second conclusion, assume that x G E>(A) an d 
x i (E). Then x G )E( . Since a(A\ E ) = a> {A) := {A : |A| > (A)} and g a> {A) 



wc 



know that A\e >( ^ a ~ ) ■ £>(A) — ^ ^>{A) is fnvertible. Let e > be arbitrary. Using the fact 
that x G ^>(A) we know that 



lirnsup ^|||A|^ w x||| < {A)~\ 

and thus there exists m G N such that for all n G N 

n > m =► ^Pj^i < (A)' 1 + e. (4.2) 
From Observation 2.1 and Theorem 2.1 we get 

x G >£><A>(r =► A\ e ^ {a x G >£> (A) ( =>■■■=► A\ e ^ {a x G >£7> (A) < , 



and from (2.2) we have 



|||x|||<)A|^ w (|P|^ w x|||, 
i <A> x|||<)Ab, w (|P|g w x|||, 



iKJx|||<>^,<lll^ (A> x 



Hence 



ll|x|||<(>A|^ (A) <n||A|^ w x|||. (4.3) 
Finally from Observation 2.1 and (4.2), (4.3) we obtain 



^ >W V IPli", .x||| (A>-i+e l + e-<A>' 



)^>>A|^ w ( r = V()^> M) <) W > f/ m.-n ..,„ > I \ \ — \ I — 



which get a contradiction with the fact that A is r-dominating. □ 

By I we denote interval [—1, lj. For e>0we put Bc(0,e) := {z G C : |z| < e}. 

Corollary 4.1. Let e > and iVeN. Assume that an operator A G £> r (K x K w_1 ) for 
r = 1/e is given. 

i) K = R. Then there exist unique eigenvalue A of A such that |A| < )A( r and the 
eigenspace is one-dimensional space. The unique (module rescaling) eigenvector x 
corresponding to the eigenvalue A satisfies 

xG (1,0,..., Of +£-(0,1,..., If. 
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ii) K = C. Then there exist unique eigenvalue A of A such that |A| < )A( r and the 
eigenspace is one-dimensional space. The unique (module rescaling) eigenvector x 
corresponding to the eigenvalue A satisfies 



x G (1,0,..., 0)' J +{0}xfi c (0, e)^ 1 C (1, 0, . . . , 0) T +e-(0, 1, . . . , I) T +e-(0, 1, . . . , ifi 

Proof. It is a direct consequence of Theorem 4.1 and Proposition 2.1. □ 
Example 4.1. Consider the following square matrix 



A 



1 0.4 0.5 
0.4 4 0.4 
0.5 0.4 8 



We find the localization of xi eigenvector corresponding to the eigenvalue of smallest 
module. Since A satisfies for r = 1.1 

)A( r < 1.81(81) < 3.24 < (A) r , 

we can use Theorem 4.1. We get 

xi G {(xi,x 2 ,x 3 ) : max{\x 2 \,\x 3 \} < —\xi\}. (4.4) 
Finding the roots of the characteristic polynomial of A, which is equal to 
W{\) = -A 3 + 13A 2 - 43.43A + 29.72 



_ (A _ 4)(A _9-v3m )(A _ £± ^m. 



we get 



/Ai-7.5 0.8Ai-3.2 \ r 

Xl= VA7^ X ' (X 1 -A)(X 1 -0.5f tX ) fOT XGM ' 

xi « (-15.686641x, 1.9070447s, x) for 
Note that xi satisfies (4.4). 

5 Applications 

In this section we show the numerical applications of our theory. For the purposes of 
this article, we use the programming language C++ with 'Boost' libraries [1] (version 
1.51.0), which contains the interval arithmetic, and software library for numerical linear 
algebra 'Linear Algebra PACKage' [2] (revision 1348). Alternatively, one could use 
numerical computing environment 'Matlab' with an interval package. In the first section 
we presented a simple algorithm (see Algorithm 1.1), which strictly localizes the position 
of particular eigenvalues and eigenvectors of a matrix. 

Now we discuss in more details the steps from the Algorithm 1.1. First we determine 
the numerical approximation of eigenvectors xi , . . . , x^ and corresponding eigenvalues 
Ai, ... j Ajv of the matrix A : C N — > C . In the step 1) of this algorithm we create matrix 
P of numerical approximate eigenvectors Xi , . . . , Xn . For the step 2) we calculate the 
inverse matrix of the interval matrix which is composed of the approximate eigenvectors. 
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The inverse matrix we can determine by the method of Gauss- Jordan elimination. In 
the next step we calculate the interval matrix 



J=[P- 1 ]-[A]-[P]. 

The most important stage of our method is the step 4). In this step we want to take 
advantage of Corollary 4.1, but first we need to check whether our matrix J satisfies 
the assumptions of this corollary (is it r-dominating) . We want to localize the fc-th 
eigenvector of a general matrix. We create matrices B £ (see (1.1)) and we check whether 
the following inequality is satisfied (see Theorem 3.1): 

}B £ (<(B £ ). 

If we find an e which satisfies the above condition for the matrix B e , then from Corollary 
4.1 we know that the k-th eigenvector of J satisfies 

(l,0,...,0) T + e-(0,I,...,I) T + e-(0,I,...,I) T i, 

where I = [—1,1]. Since J = [P -1 ]] • {A} ■ [P] we obtain the A;-th eigenvector of A 
multiplying the k-th. eigenvector of J by interval matrix [P]. Finally we have 



x k G pc k J + l-e, e_ 



N 

■E 



N 



[Re(xi)] + [-e,e].£pm(xi)]i. 

i^k 



The eigenvalue corresponding to the eigenvector x k we estimate from 

A ^ ^ (Ax k ,x k ) 
(x k ,x k ) 

Now we show some examples where we used our method. 



(5.1) 



Example 5.1. Consider 500 random matrices of size 5x5 with randomly generated 
values from the set [—1,1]. For each such matrix we launched our program to find 
accurate e location of all eigenvectors. We obtained e € [10~ 15 , 10~ 5 ] and for the results 
we have the fist-quantile to equal 6 • 10~ 13 , the median to equal 10 -11 and the third- 
quantile to equal 8 • 10~ 10 . 

Example 5.2. Let A be given in the matrix form 



A 



-2 





3 


1 


-1 


3 


-3 


-2 


4 


1 





3 


-1 


4 


2 


2 


-3 





-4 











-3 


-4 


2 


-2 





-3 


-1 


1 


-4 


-1 


4 


1 


-1 


2 


4 


1 


2 





3 


-1 


4 





4 


3 


-2 





1 


3 


4 


-1 


1 


2 


1 


-4 


2 


-2 


-4 


-2 





4 


1 


-1 


-2 


-4 


-2 


4 


1 


-1 


-3 


-4 


2 

















4 





4 


-4 


2 





-1 





-2 


-4 


4 





-4 


-3 


4 


4 





4 


-3 


3 


3 


-1 


the followin 


g set of e £ [10 


-10 


,KT 8 


]■ 
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Due to the large size of eigenvectors of A we write only first of them 



xi €( 0.0032301+ [32, 51] 
0.0032301+ [20, 62] 
0.1312400+ [61, 85] 
0.2755005+ [42, 60] 
0.1849394+ [66, 98] 
0.1849394+ [66, 98] 
0.3754802+ [56, 77] 
0.3754802+ [56, 77] 
0.1456179+ [89, 94] 



10" 



(-0.1909691 + [-14, -10] • 10~ 9 ) 



10~ 9 + (0.1909691 + [ 9, 16] • 10" 9 )i, 
10~ 9 + [-34, 34] • 10 



ht 9 + 



28, 28] • 10" 
10" 9 + (-0.0068575 + [-54, -45] • 10~ 9 )i, 
10~ 9 + (0.0068575 + [ 45, 54] • 10~ 9 )i, 
10~ 9 + (0.0266195 + [ 36, 44] • 10~ 9 )i, 
10" 9 + (-0.0266195 + [-43, -36] • 10~ 9 )i, 
10~ 9 + [-54,54] • 



-0.0358730+ [-7,0] • 10~ 9 + [-78,78] • 10" n i), 
The eigenvalue corresponding to the eigenvector xi satisfies 

Ai G 5.56625 + [46, 81] • 10~ 7 + (3.1629 + [69, 72] • 10~ 6 )i 
Example 5.3. Consider the polynomial 

W(x) = x 5 + (5 - i)x 4 - 7ix 2 + (2 + 4z)x - 8. 
To find all the roots of polynomial W(x) we estimate the eigenvalues of the matrix 



A 



From our program we obtain 







1 






1 



1 





-2 -Ai 7i 

-9 







1 

-5 + i 



Ai G -5.1189735 + [-66, -07] • 10~ 9 + 1.2610393 + [39, 67] • W~ 9 i, 
A 2 G 0.384 + [ 37, 42] • 10~ 5 + 1.2215 + [ 45, 97] • 10~ 6 i, 
A3 G 0.9572 + [ 53, 90] • lO" 6 + 0.1374 + [ 16, 42] • 10~ 6 i, 

A 4 G -1.0805 + [-80, -10] • 10~ 6 - 0.6647 + [-95, -42] • 10~ 6 i, 

A 5 G -0.1421 + [-68, -30] • 10" 6 - 0.9552 + [-97, -45] • 10~ 6 i. 
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